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^ : ABSTRACT 

o 

' The "Nessie" Nebula is a filamentary infrared dark cloud (IRDC) with a large 

aspect ratio of over 150:1 (1.5°x 0.01°, or 80 pc x 0.5 pc at a kinematic distance 
K^ \ of 3.1 kpc). Maps of HNC (1-0) emission, a tracer of dense molecular gas, made 

;-i I with the Australia Telescope National Facility Mopra telescope, show an excellent 

morphological match to the mid-IR extinction. Moreover, because the molecular 
line emission from the entire nebula has the same radial velocity to within ±3.4 
kms~^, the nebula is a single, coherent cloud and not the chance alignment of 
multiple unrelated clouds along the line of sight. 

The Nessie Nebula contains a number of compact, dense molecular cores which 
have a characteristic projected spacing of ~ 4.5 pc along the filament. The theory 
of gravitationally bound gaseous cylinders predicts the existence of such cores, 
which, due to the "sausage" or "varicose" fluid instability, fragment from the 
cylinder at a characteristic length scale. If turbulent pressure dominates over 
thermal pressure in Nessie, then the observed core spacing matches theoretical 
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predictions. We speculate that the formation of high-mass stars and massive 
star clusters arises from the fragmentation of filamentary IRDCs caused by the 
"sausage" fiuid instability that leads to the formation of massive, dense molecular 
cores. The filamentary molecular gas clouds often found near high-mass star- 
forming regions (e.g., Orion, NGC 6334, etc.) may represent a later stage of 
IRDC evolution. 

Subject headings: ISM: clouds - stars: formation 



Introduction 



When infrared dark clouds (IRDCs) were originally ide ntified as dark extinc t ion features 



against the Galactic background at mid-IR wavele ngths (jPerault et al.l Il996t lEgan et al. 



19981 : ICarey et al.lll998l . l2000l : iHennebelle et al.ll200ll ). it was not clear what role, if any, they 



might play in the process of star formation. Their high column densities (~ 10-10 cm~ ), 
together with their cold temperatures (< 25 K) suggest that if stars do form within IRDCs, 
then IRDCs must comprise an early phase in the process; otherwise the embedded stars 
would have heated their surrounding gas and dust. Recently, compact molec ular cores within 



IRDC s have been found to contain embedded young stars or protostars (e.g.. [Chambers et al. 



20091). A few of these embedded young ste l lar objects will evolve into high-mass stars (e.g. 



Beuther et al 



2005 



Rathborne et al.ll2005l . I2006L 120071 ISimon et al.ll2006l : IPillai et al.l 12006 



Beuther fc Steinackerl 120071 ). This star- formation activity, along with the similarity in size 
and mass between IRDCs and cluster-forming molecular clumps, has led to the suggestion 
that IRDCs are the b irthplaces of all high- mass stars and clusters (e.g., 
20061 : IWvrowskilboOsh . 



Rathborne et al. 



One striking feature of IRDCs is their filamentary shape, evide nt in the images from 
mid-IR Galactic surveys such as M6'X( ]Price et al.l l200lll. GLIMPSE feeniamin et al.ll2003l ). 
and MIPSGAL (jCarey et al.l 120051 ). Similar filamentary structures are also found in famil- 
iar examples of a s sociated high-mass star-forming rnolecu lar clouds, such as Orion (e.g.. 
Plume et al.ll2000t iTatematsu et al.ll2008l : iBally et al.lll987l ). Because these warmer clouds 
bear a striking resemblance to IRDCs, the warm, star-forming filamentary clouds may rep- 
resent a later evolutionary stage in the life of an IRDC. Indeed, filamentary clouds may play 
an important rol e in high-mass star formation because they allow for enhanced accretion 
rates onto cores (JBanerjee fc Pudritzll2008l : lMyersll2009l ). 



An extreme example of a filamentary IRDC is the "Nessie" Nebula centered near (/, 6) = 
(338.4°, —0.4°). Nessie is identified as a dark mid-IR extinction feature in the GLIMPSE 
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and MIPSGAL surveys (Figured]). If the mid-IR extinction represents a single object, then 
this IRDC has an unusually large aspect ratio of over 150:1 (1.5° x 0.01°). Nessie contains a 
number of extinction maxima which in other IRDCs are t ypically associated with pre-stellar 



and protostellar cores with masses of order 100 Mq (e.g., iRathborne et al.ll2006l ) 



If the formation of IRDCs, perhaps due to the passage of a spiral shock, results in fila- 
ments, then these long, approximately cylindrical structures may produce the necessary con- 
ditions for the formation of dense molecular cores, and ultimately of clusters and high-mass 
stars. Indeed, the theory of gravitationally bound gaseous cylinders predicts the formation 
of such dense cores, which due to the "sausag e" or "varicose" fluid instability, fragment frorn . 



the cylinder at roughly regular intervals (e.g., Chandrasekhar &: Ferm|l953l : lNagasawalll987l : 



Inutsuka fc Miyama]ll992l : lNakamura et al 



1993 



Tomisaka 



1995 



19961 ). This instability pro- 



vides an attractive physical mechanism for high-mass star formation and cluster formation 
in filamentary IRDCs. 

In this Letter we report on molecular line mapping of the Nessie Nebula. We compare 
our results with the theory of self-gravitating fluid cylinders and show that the observations 
broadly conform to the expectations of the "sausage" instability. We speculate that this 
mechanism may be important in the formation of high-mass stars and star clusters. 



2. Observations 



Millimeter molecular line observations of the Nessie Nebula were taken on 2008 August 
with the Australia Telescope National Facility (ATNF) 22 m Mopra Telescope. The backend 
consisted of the wideband MOPS correlator configured to observe simultaneously 16 sepa- 
rate passbands covering the frequency range 86 to 93 GHz. Each of these passbands was 
137.5 MHz wide, with 4096 spectral channels, corresponding to a velocity resolution of 0.11 
kms~^. We chose individual passbands to include molecular lines of interest. The brightest 
of the detected lines was HNC (1-0) (90.66 GHz). Although we detected a number of fainter 
molecular lines, these yield similar results and will be presented in a subsequent paper. 

The nebula was mapped in the "on-the-fly" mode, in which data are collected while 
the telescope scans in a raster pattern. The raster rows were chosen to be perpendicular 
to the Galactic plane in order to avoid striping artifacts parallel to the nebula's long axis. 
The spacing betw een rows was 12", or about 1/3 of the Mopra beamwidth of 36" FWHM 
(JLadd et al.ll2005l ). Flux calibration with a noise diode and an ambient temperature load 
was performed every 20 minutes. 



Data reduction and analysis utilized the ATNF livedata and gridzilla programs. The 
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individual spectra were co-added and gridded onto a uniform 15" grid. Linear baselines were 
removed from line-free channels. All spectra are presented on the antenna temperature (T^) 
scale. To convert to main-beam brightness temp eratures, orie sho uld divide the antenna 
temperatures by the main-beam efficiency of 0.49 (JLadd et al.ll2005l ). System temperatures 
for the observations were ~ 180 — 300 K, which yielded a typical rms noise of T\ = 0.14 K 
in each spectral channel. 



3. Results 



The Mopra integrated intensity HNC (1-0) maps closely correspond to the regions of 
mid-IR extinction (Figured]). The map consists of fainter, uniform emission associated with 
the filament, upon which compact, brighter regions we call "cores" are superposed. Gaussian 
fits to the HNC (1-0) lines demonstrate that every position within Nessie has essentially the 
same radial velocity, —38 ±3.4 kms~^ (Figured]). If the Clemens rotatio n curve is used, 
this velocity corresponds to a kinematic distance of 3.1 kpc ( 1Clemenslll985l ). 



The molecular cores in Nessie are often associated star-formation activity. Specifically, 
many cores show the p resence of excess 4.5 Mm ern i ssion called "green f uzzies," which indi- 
cates shocked gas (e.g., iNoriega-Crespo et al.ll2004l : iMarston et al.ll2004J). and bri ght 24 /xm 



point sources, which indicate an embedded protostar (e.g.. [Chambers et al.ll2009l ). 



The cores do not appear to be randomly spaced within the filament, but instead have a 
characteristic spacing. The filament is thus reminiscent of a string with beads that are spaced 
sparsely but approximately uniformly along its length. We used the c lumpfind2d algorithrn 
with a threshold of seven times the rms noise and 2.5a increments (j Williams et al.l Il994j ) 
on the HNC (1-0) integrated intensity map to identify 12 molecular cores (see Figured])- 
The exact choice of clumpfind parameters will result in the identification of slightly more or 
fewer cores, but the resulting core spacing is quite insensitive to reasonable input parameters. 
Since the cores are relatively isolated and well separated, clumpfind works well and identifies 
the cores that one would tend to select by eye. The positions and selected properties of the 
cores are listed in Table dl Columns 1 and 2 give the Galactic coordinates of the position of 
peak integrated HNC intensity in each core. Column 3 the peak HNC integrated intensity. 
Columns 4 and 5 the angular size and physical extent. Column 6 the LTE mass (see Section 
H]), and Column 7 whether or not the core contains a 24 /im point source. The mean projected 
spacing between cores is roughly 4.5 pc. 

Although we prefer to use the molecular gas clumps as our indicator of core spacing, 
since we thereby include pre-stellar cores, another useful indicator of cluster-forming sites 



-5- 



would be 24 /iin point sources, which indicate embedded high-mass stars or protostars. 
Indeed, about half of the molecular cores contain 24 /im point sources (see Tabled]). For all 
of the unresolved 24 /im sources associated with the nebula, we find typical spacing of 4.1 
pc, a value within 10% of the molecular core spacing. Thus, both the spacing of gas clumps 
and of 24 /im sources are in satisfactory agreement. 



Discussion 



The existence of cores that are spaced at approximately a characteristic spacing within 
a highly filamentary cloud is in broad agreement with the predictions of fragmentation of a 
self- gravi tating fluid cylinder due to the "sausage" instability. The theory's original presen- 
tation by IChandrasekhar fc Fermil (119531 ) treated an incompressible fluid. Later refinements 
included an isothermal, thermally supported cylinder with various magnetic field configu 



ratio ns (JNagasawal 119871 : llnutsuka fc Miyamal Il992l : iNakamura et al.l Il993l : iTomisakal Il995 
1996h . 



Although the fragmentation due to the "sausage" instability of a self-gravitating fluid 
cylinder is similar to three-dimensional Jeans collapse, there are important differences. In 
spherical gravitational Jeans collapse, perturbations of all wavenumbers grow at the same 
rate, but in cylindrical gravitational collapse perturbations with certain wavenumbers grow 
more quickly than others. Perturbations with the most unstable wavenumber will grow 
fastest, and hence overdensities (cores) will tend to form at a length scale whose wavenumber 
corresponds to this fastest growing mode. For a cylinder of infinite radius, this len gth scale is 
maintained even in the presence of a magnetic field parallel to the filament's axis (JNagasawa 



1987h . 



The theory thus predicts (1) that multiple cores should form within a filament, and (2) 
that the spacing between these cores should be roughly periodic, with a characteristic length 
scale equal to the wavelength of the fastest growing unstable mode of the fluid instability. 
For an incompressible fluid, th is wavelength is Xmax = H-R where R is the cylinder's radius 
(IChandrasekhar fc Fermilll953l ). In an infinite isothermal gas cylinder, the relation is Xmax = 
22H, where H is the isothermal scale height given hy H = Cs(4:TiGpc)^^^'^, where Cg is the 
sound speed, G th e gravitational constant, and Pr the gas ma ss density at the center of the 
filament {R = 0) (JNagasawal 119871 : llnutsuka fc Miyamal 1 19921 ). For isothermal cylinders of 
finite radius R surrounded by an external, uniform medium, the spacing depends on the 
ratio between the cylinder radius and the isothermal scale height, R/H. For R ^ H the 
core spacing approaches that for an infinite radius cylinder, Xmax = 2211, but ioi R <^ H 
the spacing reduces to that of the incompressible case, Xmax = H-R- 
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We now compare the observations of Nessie to this theory under the assumption that 
Nessie is well approximated by an isothermal cylinder. If Nessie lies in the limit R -^ H, the 
theory reduces to that of an incompressible fluid, and predicts cores spaced at a characteristic 
length scale ~ Hi?. One can estimate the radius R of Nessie from the radial extent of the 
sharp extinction edges evident in the GLIMPSE and MIPSGAL im ages. This su ggests a 



radius of i? ~ 0.01° or ~ 0.5 pc at the kinematic distance of 3.1 kpc (jClemenslll985l ). Thus, 
one would expect cores to form at a spacing of ~ 5.5 pc, in satisfactory agreement with the 
observations if the filament axis is mostly perpendicular to the line of sight. 

If Nessie is treated as an isothermal gas cylinder with a central volume density n = 10^ 
cm~^ and T = 10 K, appropriate for cold, dense, molecular gas traced by HNC (1-0) emission, 
then the isothermal scale height is if = 0.05 pc. Thus, Nessie would be in the regime R^ H, 
and the theoretical spacing between cores should be Xmax = 22H ~ 1 pc. This is smaller 
than the observed spacing by about a factor of 5. This discrepancy may arise from the fact 
that turbulent pressure dominates over thermal pressure. 

Theoretical studies of self-gravitating cylinders usually assume that thermal pressure 
is the dominant gas pressure. However, because the observed linewidths in most molecular 
clouds typically exceed the thermal sound speeds by large factors, turbulent pressure usu- 
ally dominates over thermal pressure. If one rep laces the sound speed Cg with the velocity 



dispersion a to account for turbulent pressure (iFiege fc Pudritzll2000l ). the effective scale 
height Heff is increased considerably over the thermal scale height H. Assuming a central 
volume density again of n ~ 10^ cm~^, and a FWHM linewidth of 2.5 kms~^, typically ob- 
served toward Nessie, the effective scale height is ife// ~ 0.17 pc, which leads to a spacing of 
Xmax = '^'^Heff ~ 4 pc between the cores. Again, the theoretical prediction is in satisfactory 
agreement with the observations, especially since the central density and hence the scale 
height H is uncertain. Moreover, the theoretical assumption of a uniform, isothermal cylin- 
der (sometimes embedded in a uniform pressure external medium) is obviously idealized, 
and the fragmentation length scale may well differ in more realistic treatments. 

The theory also predicts that self-gravitating cylinders in equilibrium have a maximum, 
critical linear mass density (mass per unit length along the cylinder's axis). Above the critical 
value, the cylinder would collapse radially into a line. This critical linear mass density is 
given by {M/l)max = 2v'^/G; where v is the sound speed Cs in the case of thermal support 
(Stodolkiewicz 196 3, Ostriker 1964) and t he turbulent velocity dispersion a in the case of 



turbulent support ( jFiege fc Pudrita 120001 ). If turbulent pressure dominates over thermal 
pressure in Nessie, the critical linear mass density can be estimated as {M/l)max = 84(AV^)^ 
Mq pc~^, where AV^ is in units of kms~^. For typical line widths observed in Nessie {AV ~ 
2.5 kms~^), {M/l)max ~ 525 M© pc^"*^. [Note that magnetic fields can alter the critical 
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linear mass density, increasing it for poloidal fields an d decreasing it for toroidal fields, but 
probably by small factors ~ 1 f lFiege fc Pudrita 120001 ) .] 



The mass per unit length in the Nessie Nebula can be estimated from the HNC (1- 
0) fluxes, although this estimate is difficult due to uncertainties in molecular abundances, 
depletion of molecules frozen onto cold dust grains, and optical depth effects. Nevertheless, 
assuming optically thin HNC (1-0) e mission, LTE, a tena perature of 10 K, and an HNC 
abundance of y^HNC/H2 = 2 x 10~^ (JTennekes et al.l l2006l ) appropriate for cold clouds, we 
arrive at a typical linear mass density in Nessie of 110 Mq pc~^. Because this value is smaller 
than the maximum predicted by theory, the estimated linear mass density is consistent with 
theoretical expectations. Unfortunately, the reported abundan ces of HNC vary wid ely, with 
values for cold clouds ranging from 4.4 x 10^^^ to 2 x 10~^ (jWootten et al.l 119781 ). These 
HNC abundance values give rise to a wide range in the linear mass density of up to 5 x 10^ 
Mq pc~^. Future observations of the dust continuum will better constrain the linear mass 
density. 

It is interesting to compare the value of M/l obtained by theory and estimated in Nessie 
with that measured for the molecular fila ment associated w ith Orion, M/l = 385 Mq pc^^, 
assuming M = 5 x 10'^ Mq and / = 13 pc (JBally et al.lll987l ). This value is within a factor of 
a few of that estimated for Nessie, and also smaller than the maximum value predicted by 
theory. 

The maximum mass of the cores can also be estimated from theory. The core masses 
should roughly be given by M ~ Xmax{M/l). Assuming the maximum, critical value for 
M/l of ~ 525 Mq pc"^ for the case of turbulent pressure support, and a spacing of 4.5 pc, 
one expects the cores to have maximum masses Mmax ~ 24OOM0. If instead the average 
linear mass density of 110 Mq pc~^ derived for Nessie is used, the core masses should be no 
larger than about 500Mq. For Nessie, the estimated core masses based on the integrated 
HNC (1-0) fiuxes over the solid angles range from M ~ 30 - 39OM0 (Table [I]). 

To summarize, if turbulent pressure, rather than thermal pressure, provides the support 
for Nessie, the core spacing (4.5 pc) agrees with theoretical predictions. Moreover, both the 
linear mass density and the core masses in Nessie are also consistent with the theory, since 
they are comparable to, but fall below, the maximum theoretical values. 

We speculate that all cluster formation arises from the fragmentation of filamentary 
IRDCs due to fiuid instabilities that lead to massive, star-forming cores. The periodic spacing 
of such cores within Nessie strengthens this suggestion. Moreover, filamentary molecular gas 
clouds are usually associated with high-mass star-forming regions, such as Orion, and these 
warmer clouds bear a striking resemblance to IRDCs. Such clouds may represent a later 



evolutionary stage of a filamentary IRDC 

A cartoon sketch of the proposed evolution of IRDCs is given in Figure [H Initially a 
filamentary IRDC is formed. The filament reaches an equilibrium configuration with a linear 
mass density close to the critical value. The filament then fragments due to the "sausage" 
instability into roughly periodically spaced cores. Such dense molecular cores have sufficiently 
large masses to form clusters. The mass of the cores is set by a combination of the turbulent 
velocity dispersion, which determines the linear mass density, and the wavelength of the 
fastest growing mode of the "sausage" instability. Thus, the cluster-forming core mass, and 
ultimately, the mass of the resulting cluster, is determined by the physics of the turbulent, 
self-gravitating cylindrical fiuid. Within the more massive cores, high-mass protostars form. 
When an embedded O or an early B protostar eventually reaches the main sequence, it ionize s 



the surrounding gas and creates a bubble in the interstellar medium ( IChurchwell et al.ll2006l ) . 
Eventually the embedded cluster emerges. 

In the earliest phases, then, filamentary IRDCs with several cores would be observed. 
As the H II regions form, bubbles would be associated with the IRDCs, and eventually the 
IRDC would be segmented and disrupted as the Hii regions grow and overlap. Later, the 
filament would contain one or more bright clusters, and appear similar to the filamentary 
molec ular clouds associate d with high-mass star-forming regions such as Orion and NGC 



6334 (JKraemer et al.lll999[ ). Examples of each of these proposed stages can easi ly be found, 



from dark filamentary IRDCs like Gil. 11 — 0.11, sometimes called the Snake (jCarey et al. 



19981 . I2OOOI ). to IRDCs containing one or more bubble-like H 11 regions such as Nessie, to 



active high-mass star-forming regions like Orion and NGC 6334. 



5. Conclusions 

We have imaged a filamentary IRDC, the Nessie Nebula, in HNC (1-0) molecular line 
emission with the ATNF Mopra Telescope. Because the HNC (1-0) radial velocities are 
uniform across the filament, Nessie is a coherent, single cloud. Both HNC (1-0) and mid-IR 
images reveal that Nessie contains numerous dense cores, which appear to have a roughly 
periodic spacing of 4.5 pc. 

The theory of self-gravitating, gaseous cylinders shows that they are unstable to a fiuid 
instability, called the "sausage" instability, which causes the filament to fragment into cores 
with a roughly constant spacing. Thus, unlike the Jeans collapse of uniform media, the 
fragmentation of cylinders has a preferred length scale. Theory predicts that the fastest 
growing wavelength scales as the radius R or the isothermal scale height H of the cylinder. 
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depending on the ratio R/H. In the conditions thought to be appropriate to Nessie, the 
spacing between clumps is predicted to be 4 to 6 pc, in good agreement with the observed 
spacing of 4.5 pc. 

Because of the ubiquity of the association between filamentary molecular clouds and 
high-mass star forming regions, it is tempting to suggest that the physics of molecular 
filaments has a profound influence on cluster formation. Indeed, such filaments may be a 
necessary initial condition to form dense cluster-forming cores with the observed properties. 
We speculate that the "sausage" instability may be the dominant physical mechanism to 
produce cores from filamentary IRDCs. We suggest an evolutionary sequence in which 
filamentary IRDCs fragment into cores, which in turn spawn star clusters and high-mass 
stars. Thus, familiar high-mass star-forming clouds such as Orion may have begun their 
lives as filamentary IRDCs. 

The authors gratefully acknowledge the funding support through NASA grant NAG5- 
10808 and NSF grants AST-0098562, AST-0507657, and AST-0808001. The authors are 
grateful to Prof. Mark Johnson of Northwestern University for suggesting the idea of varicose 
instabilities, to Julia Duval-Roman for her assistance, and to an anonymous referee for 
important suggestions. 
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Fig. 1. — Top panel: a false three-color image of the Nessie Nebula. The 3.6 /im (blue) and 
8.0 /im (green) emission is from GLIMPSE, and the 24 ytim (red) emission is from MIPSGAL. 
Second panel: the false three-color mid-IR image of the Nessie Nebula from the top panel 
overlaid with integrated HNC (1-0) contours from the Mopra telescope. Note the excellent 
correspondence between the HNC emission and the 8 /zm extinction. Third panel: an HNC 
(1-0) integrated intensity map from the Mopra telescope with core positions marked with 
cyan crosses. Bottom panel: a velocity-field (first moment) map from the HNC (1-0) map 
of the Nessie Nebula. All of the molecular emission is at the same velocity to within ±3.4 
kms~^, demonstrating that the filament is a single coherent object. 
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Fig. 2. — Hypothetical sequence of IRDC evolution. Time increases downward in the figure 
from top to bottom showing different evolutionary stages: a filamentary IRDC is formed, the 
filament fragments into dense cores, massive stars and Hii regions form within these cores, 
and finally star clusters emerge. 
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Table 1. Properties of Cores in the Nessie Nebula 



Coordinates 


Peak Integrated HNC T^ 


Angular Size 


Physical Extent 


Mass 


24 /im Point Source? 


1° 


6° 


(K km s- 


■i; 


( arcmin ) 


(pc^) 


(Mq) 




338.17 


-0.47 


1.15 




3 


2 


80 


Y 


338.32 


-0.41 


2.19 




7 


6 


250 


Y 


338.41 


-0.41 


2.26 




9 


7 


320 


Y 


338.46 


-0.43 


1.03 




2 


1 


50 




338.63 


-0.46 


1.42 




4 


4 


130 




338.68 


-0.46 


0.99 




1 


1 


30 




338.77 


-0.46 


1.14 




2 


1 


40 




338.8 


-0.47 


1.2 




4 


4 


110 




338.87 


-0.48 


2.24 




6 


5 


220 




338.94 


-0.49 


3.05 




10 


8 


390 


Y 


339.02 


-0.41 


1.47 




7 


6 


220 




339.06 


-0.41 


2.02 




9 


7 


310 


Y 



